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^ \ We discuss the impact of different measurements of the d/u asymmetry in 

the extraction of parameterizations of parton distribution functions. 
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I. INTRODUCTION 



Recently, the E866/NuSea Collaboration |T[ has released the final results corresponding 
to the analysis of their full data set of Drell-Yan yields from an 800 GeV/c proton beam 
on hydrogen and deuterium targets. This analysis, aimed to size the d/u asymmetry in the 
proton, has confirmed previous measurements by E866 0, but with a substantial increase 
in precision and extending the kinematical coverage. 

The measurement of the d/u asymmetry in the proton has been pursued in recent years 
by several groups , triggered by the NMC Collaboration measurement of deep inelastic 
muon scattering from hydrogen and deuterium targets, which suggested that the d quark 
^Cj . content in the proton was larger than that of the u quarks ||. 

The ongoing interest in this issue has been not only the verification a hitherto unexpected 
flavour symmetry breakdown, presumably of non-perturbative origin, but also the precise 
determination of the ratio d/u, as required by modern global QCD fits of parton distribution 
functions [PHTllI- These global analyses are an essential ingredient in the perturbative QCD 
description of hadrons and their hard interactions, and in recent years have evolved in 
accuracy and sophistication along the increase in precision of data |HZ|. 



Indeed, prior to NMC finding, global extractions of parton distributions were made 
under the assumption of flavour symmetry for the light sea quarks (u = d). Most recent 
sets, however, were constrained to satisfy either NMC or, more likely, E866 data, imposition 
that also influences the extraction of valence quark densities P-pH. 

It is interesting to notice at this point, that the constraints on sea quark densities coming 
from the various asymmetry measurements are not necessarily equivalent and, eventually, 
they influence in different ways the extraction of parton densities in global fits. These 
differences are most noticeable when comparing for example E866 Drell Yan data and NMC 
results on the difference between the proton and the neutron structure functions. In other 
words, global fits designed to accommodate E866 data fail to reproduce that coming from 
NMC, and conversely. 
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Given the accuracy of modern parton distributions, and the fact that F$ arid _F 2 d , as well 
as E866 data are nicely accounted for by them, it is somewhat surprising that they fail to 
reproduce F$ — data at intermediate values of the variable x, as shown in Figure 1. 
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FIGURE 1: NMC (a), E886 (b), and HERMES (c) d/u- related measurements against current 
LO parameterizations. 

The failure in reproducing the different data sets simultaneously may be pointing out the 
need of more precise data, the necessity of a revision of some of the underlying assumptions 
in the analysis of the data, such as the role of nuclear corrections in the use of deuterium 
targets flT3|-(r5|], for example, or even the presence of some other hitherto neglected effect. 
In the mean time, it is worthwhile analysing to what extent parton distributions are affected 
by the different alternatives. 
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In order to assess the impact of different results on the d/u asymmetry in the extraction 
of parton distribution functions, we perform in this paper a leading order (LO) QCD global 
extraction of parton densities from a large set of DIS data as usual, but forcing the agreement 
with one particular flavour asymmetry measurement (NMC, E866, or HERMES) as an ad 
hoc constraint in the fit. This is accomplished giving to the 'preferred' data set a larger 
statistical weight than it actually has, and it is done also in order to increase the little 
relative statistical relevance that they have against the bulk of DIS data. 

As result of this approach we find that even though one can produce parton densities 
using either E866 or NMC data to implement the flavour symmetry breaking, the constraints 
that these sets of data impose on the distributions can not be satisfied simultaneously in 
any way. HERMES data show a similar conflict with NMC, but is in rather good agreement 
with E866. 

Parton distributions in agreement with E866/HERMES symmetry breaking scenario lead 
in the fits to lower global x 2 values than those that reproduce NMC data, and the shape 
of the sea and valence quark densities of these extreme scenarios show clearly noticeable 
differences. 

In the following section we summarize the main features of the parton distributions pa- 
rameterizations used, the Q 2 -evolution strategy, and the data included in the global analyses. 
Then we present the outcome of the three fits in which either E866, NMC, or HERMES 
data have been favored, against the standard one in which no additional weights have been 
applied. Finally, we present our conclusions. 

II. PARTON DENSITIES AND DATA. 

In order to parameterize quark densities at the initial scale Qq we adopt the standard 
functional dependence as in fllO"||: 



xu v {x,Ql) = N u x au (l-x)^(l+ luV ^+5 u x) 

xd v (x,Q 2 ) =N d x ad {l-x) /3d (l + ldV ^ + 5 d x) (1) 

for valence quarks; 

x£(x,Q2) = jV E x as (l -x) fe (l + 7 SV ^ + 5 E x) (2) 
for the sum of light sea quarks and antiquarks {u s + u s + d s + d s + s s + s s ) and 

xg(x,Q 2 ) = NgX as (l-x)^(l+7 g y/x + 5 g x) (3) 



for gluons. In order to be consistent with CCFR data on dimuon production fl6fl , strange 
sea quarks are assumend to be s = 0.1 £ while u and d densities are given by 

x(d + u) = 0.4 £ 

x(d-u) = N A x aA (l - xf A (l + -f A ^ + 5 A x 2 ) (4) 

Leaving aside the normalizations for valence quarks and gluons, which are fixed by charge 
and momentum conservation, the above distributions imply 21 free parameters to fit. 
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The initial scale Qq is taken to be 1 GeV 2 and the evolution is performed using a varia- 
ble number of active flavours with thresholds at 1.69 GeV 2 for charm, and 18.49 GeV 2 for 
bottom. Bellow these thresholds the contributions to the structure function are given by the 
corresponding photon-gluon fusion diagrams. For the strong coupling a s the LO expression 
with Aqcd = 0.174 GeV (for four flavours) is taken. 

The data set used in the global fitting procedure is listed in Table 1, where we have 
excluded data with x > 0.75, x < 0.001, Q 2 < 2 GeV 2 , and W 2 < 7 GeV 2 in order to avoid 
higher order QCD and mass effects. 



Process 


Experiment [Ref.] 


Observable 


#Data 


DIS 


BCDMS [f| 




(177,159) 


£665 |J 


F$,Fi 


(53,53) 


HI \)^,ZEUS @ 


F$ 


(150,158) 


NMC (§ 


F$,Fj 


(130,130) 


NMC @ 


Ff - £ 2 " 


(12) 


SIDIS 


HERMES § 


d—u 
it — d 


(5) 


Drell-Yan 


£866 § 


<?pd 
2 (7pp 


(15) 


W + W~ Asymmetry 


CDF [fl| 




(11) 



TABLE I. Sets of data. 



Notice that in addition to NMC data on the proton and the deuteron structure functions, 
F% and F$ respectively, we have also included the data on the difference — F% as reported 
in reference [[J , which we take as an alternative indicator of the flavour asymmetry, in the 
same footing as E866 or HERMES data. 

III. SCENARIOS 

As it was stated above, it is clear that the last four sets of data in Table 1 have very 
little statistical weight relative to DIS data. In order to preserve the information contained 
in them, and highlight the alternative scenarios that they imply, we give them an artificial 
extra weight in the x 2 minimization procedure. For E866, NMC data on the difference 
F% — F£ and HERMES, these weights are varied according to which experiment we want 
to particularize, yielding three distinctive parton distribution sets labeled Set I, Set II, and 
Set III, respectively. For comparison we also include a set in which no additional weights 
are given, labeled as Set 0. 

In Set I, a moderate weight factor (cue866 = 10) is applied to E866 data, allowing an 
excellent matching (xi;866 = 13.4) shown by the solid line in Figure 2b., a rather good 
agreement with HERMES (Xhermes = 5-3 Figure 2c), while a rather poor accord with 
NMC (x 2 nmc = 1^3 Figure 2a.). This set has a close resemblance with MRST, and clearly 
favors the agreement with E866 data relative to that with NMC. However, at variance with 
MRST, because it includes explicitely NMC data, it also reproduces it slightly better. It is 
also very similar to the unweighted fit, labeled as Fit 0, due to the fact that E866 data has 
by itself the largest constraining power of the three sets of data. In Table 2 we show the 
corresponding x 2 values obtained. 
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FIGURE 2: NMC (a), E886 (b), and HERMES (c) d/u-related measurements against fits. 

Set II, shown as dashes in Figure 2, have also a moderate weight factor on MNC data 
(ujnmc = 8), which improves dramatically the agreement with this data (xnmc = 45), 
however E866 is not reproduced at all (xesgg = 142). The agreement with HERMES is 
also significantly reduced (Xhermes = 16)- Notice also that the overall quality of the fit 
deteriorates in a 10 %. 

Set III is designed to accommodate HERMES data {Xhermes = 4-6) with a very strong 
weight factor {ujhermes — 180), however it also reproduce fairly well E866 data (xls66 = 
19). The agreement with NMC is completely lost (Xnmc = 184), and it is even worse than 
in Set I. 

It is interesting to notice that the most striking difference between the sea quark densities 
inspired in the E866/HERMES and NMC scenarios, respectively, is that whereas NMC data 
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favors d quark densities up to more than twice larger that u ones, both E866 and HERMES 
lead to not so large d/u ratios, as shown in Figure 3a. There is also a clear difference in 
the value of momentum fraction at which the ratio reaches its maximum, which in turn 
determines the position of the peak of the d — u distribution. Both HERMES and E886 lead 
to peaks around x ~ 0.1 while NMC favors x ~ 0.2 as can be seen in Figure 3b. 




FIGURE 3: a) d/u and b) d — u distributions from the fits. 
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FIGURE 4: Valence distributions from the fits. 
Given that neutral current structure function data are only sensitive to the sum of 
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quark- ant iquark densities of each flavour, the size of the flavour symmetry breaking in the 
sea introduces an additional uncertainty in the extraction of up and down valence quark 
densities, as shown in Figures 4a and 4b, respectively. There, it can be seen that again the 
NMC symmetry breaking scenario leads to an interesting difference in the behaviour of the 
valence densities, most noticeable in the case of the down quark distribution. 





I (E866) 


II (NMC) 


III (Hermes) 


(wj) 


NMC if 


209.3 


177.6 


205.08 


201.6 


NMC F 2 D 


173.0 


200.5 


161.9 


167.9 


BCDMS Ff 


175.5 


256.7 


171.508 


179.7 


BCDMS Ff 


176.2 


222.1 


174.79 


176.3 


£665 £ 2 P 


61.9 


62.9 


61.8 


62.4 


£665 £ 2 D 


51.0 


53.3 


49.5 


49.6 


HlFg 


114.6 


118.7 


127.4 


116.5 


ZEUS F$ 


243.2 


239.5 


246.5 


239.9 


CDF A 


33.9 


12.5 


34.7 


27.9 


£866 a pD /2a pp 


13.4 


142.9 


19.2 


22.2 


NMC £ 2 ™ - Fl 


153.9 


44.9 


184.35 


136.8 


HERMES 


5.3 


16.2 


4.6 


6.5 


X 2 (unweighted) 


1411.3 


1547.8 


1442.1 


1389.6 



TABLE II. x 2 values obtained for different sets of data. 



IV. CONCLUSIONS 

We have analysed the consequences of different measurements of the ratio d/u asymme- 
try in the extraction of parton distributions. We have found that even though E866 and 
HERMES results on d/u lead to consistent extractions of parton distributions, NMC data 
on the difference of the proton and the neutron structure function can not be contained in 
this picture. 

In view of the next generation of high precision hadron collision experiments, the issue 
of uncertainties in parton distribution functions, and their precise origin, is a highly relevant 
one, that will certainly require further attention in the near future. 

Nuclear effects when using deuterium targets is a more or less obvious source of uncer- 
tainties in the flavour asymmetry measurements we have deal with. But at variance with the 
most conventional estimates for nuclear effects in deuterium, which lead to small corrections 
in the NMC observable at very small and large momentum fractions, the discrepancies be- 
tween parton densities inspired in NMC data and those in agreement with E866 or HERMES 
data are concentrated at intermediate values of x. 

These results suggest that even though the breakdown of the isospin symmetry in the 
sea of nucleons is indeed one of the reasons for the Gottfried sum rule deficit observed by 
NMC, it is far from being the answer to the question raised almost ten years ago. 
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VI. APPENDIX 





I (E866) 


II (NMC) 


III (Hermes) 


(wj = 0) 




0.305 


0.372 


0.314 


0.311 


Pu 


3.157 


3.114 


3.192 


3.176 


In 


-0.132 


-0.356 


-0.138 


-0.002 


Su 


16.84 


10.07 


16.94 


17.12 


ad 


0.171 


0.081 


0.168 


0.182 


Pd 


3.598 


3.394 


3.550 


3.639 


Id 


16.69 


878.9 


22.81 


22.38 


s d 


30.47 


291.5 


30.95 


33.65 


a g 


0.467 


0.778 


0.670 


0.583 


Pg 


20.59 


27.11 


20.69 


20.70 


la 


-3.502 


-2.130 


-4.165 


-3.467 


s 9 


8.084 


15.35 


8.577 


7.877 




0.484 


0.659 


0.421 


0.423 


as 


-0.175 


-0.179 


-0.210 


-0.210 




6.948 


6.377 


6.945 


6.966 


7s 


-1.086 


-3.806 


-1.012 


-1.066 


<5s 


11.22 


13.29 


11.45 


11.32 


iV A 


5.545 


41.85 


4.832 


5.557 


OA 


1.826 


3.885 


1.729 


1.913 


Pa 


10.83 


20.41 


10.94 


10.86 


7A 


11.74 


274.1 


11.31 


12.26 


Sa 


-43.65 


434.5 


-40.96 


41.03 



TABLE III. Parameters for the different sets. 
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